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چکیده 


روندیابی سیلاب یکی از الزامات مهم در مطالعات مهندسی رودخانه محسوب می‌شود. روندیابی هیدرولوژیکی در رودخانه‌های شریانی و 
رودخانه‌های فاقد آمار حوضه میانی متداول است. به این منظور نیاز به تهیه مقاطع عرضی و تعیین شیب‌ها در کلیه بازه‌های رودخانه می‌باشد. روش 
ماسکینگام می تواند با استفاده از آن ضمن صرفه‌جویی در زمان و هزینه, اطلاعات مربوط به عمق و دبی جریان سیلابی را در هر زمان مشخص نماید. 
کاربست روش‌های فراکاوشی نتایج رضایت بخشی را در این زمینه تاکنون نشان داده است. از این رو در این پزوهش, به ارزیابی کارایی الگوریتم 
جستجوی موجودات همزیست (805) در تخمین پارامترهای بهینه مدل غیرخطی ماسکینگام پرداخته شد. به‌منظور بررسی میزان مطلوبیت یافته‌های 
پژوهش, نتایج حاصل از الگوریتم موجودات همزیست (8009)» با نتایج سایر روش های فراکاوشی شامل الگوریتم وراثتی (0۸)) الگوریتم ازدحام ذرات 
(۳50)» الگوریتم رقابت استعماری (16۸) مقایسه گردید. در الگوریتم پیشنهادی» روش تابع جریمه غیرمستقیم در مدل برای جلوگیری از منفی شدن 
خروجی و ذخیره اعمال شده است. الگوریتم مذکور بهینه سراسری یا نزدیک سراسری را بدون در نظر گرفتن مقادیر اولیه پارامترها با همگرایی سریع 
پیدا می‌کند. نتایج الگوریتم 505 برای دو رودخانه ویلسون و کارده نشان دهنده کمینه‌سازی مجموع مربعات باقیمانده‌ها (9900) می‌باشد که برای 
رودخانه ویلسون با 11513 (۵/۸۵) و 350 (۱۲۸/۷۸) و رودخانه کارده با ۱155 (۰/۵۰۵) و 350 (۴/۵۵۲) می باشد و مانشد الگوریتم‌های ۳90 و 


۸ عملکرد بهتری نسبت به الگوریتم ۸) داشته است در نتیجه الگوریتم پیشنهادی می‌تواند با اطمینان خوبی به‌منظور برآورد مقادیر بهینه 


پاراتقدهای سل ماشکیتگام غیرطی,موزد اسفاده قران کیرد 


واژه‌های کلیدی: بهینه یابی» تابع جریمه غیرمستقیم» رودخانه شریانی» روندیابی هیدرولوژیکی, مدل غیرخطی ماسکینگام 


مقدمه 

سیل یک پدیده طبیعی است که می‌تواند خسارات جانی 9 مالی 
بی‌شماری در مناطق شهری» صنعتی و کشاورزی به‌وجود آورد. از این 
رو شناخت ویژگی‌های سیل نظیر مقدار اوج آن در زمان‌ها و مقاطع 
پیش‌بینی این ویژگی‌ها و به‌طور کلی تغیبرات عمق و دبی جریان در 
طول رودخانه‌ها با روندیابی سیل امکانپذیر می‌باشد. روندیابی سیلاب 


۸ ۰۲ ۲ و ۴- به‌ترتیب دانشجوی دکتری سازه‌های آبی. دانشجوی کارشناسی ارشد 
سازه‌های آبی. دانشیار و استاه گروه علوم و مهندسی آب» دانشکده کشاورزی» 
دانشگاه فردوسی مشهد. مشهد. ایران 
(#:- نویسنده مسئول: 6۴ص ۵) زاتمطرو۴ :[تمحصط) 


10۵]: 10.22067/[597.۷73412.85190 


پایین‌دست توسط هیدروگراف معلوم جریان بالادست تعیین می‌گردد. 
این عملیات برای مدیریت صحیح و طراحی دقیق بسیاری از 
پروژه‌های زیست محیطی و منابع آب مورد نیاز است (۲۳). روندیابی 
جریان یک تجزیه وتحلیل برای ردیابی جریان از طریق یک سیستم 
هیدرولوژیکی با توجه به ورودی است. روش‌های روندیابی سیل را 
می‌توان به‌طور کلی به دو دسته روندیابی هیدرولیکی (گسترده) و 
روندیابی هیدرولوژیکی (متمرکز) طبقه‌بندی کرد (۶). روش‌های 
روندیابی هیدرولیکی بر مبنای حل عددی معادلات سنت ونانت و 
معادلات غیردائمی متغیر تدریجی در کانال‌های باز می‌باشد. درحالی 
که روش‌های هیدرولوژیکی بر پایه حل معادله پیوستگی و روابط بین 
دبی و ذخیره موقت حجم اضافی آب در دوره سیلاب استوار است. 
روش های هیدرولوژیکی درمقایسه با روش‌های هیدرولیکی دارای 
محاسبات ساده و سریع‌تر هستند و در کاربردهای مهندسی با اطمینان 
قابل قبولی به‌کار می روند. پرکاربردترین و مشهورترین روش در این 
دسته. روش ماسکینگام است (۷). روش ماسکینگام اولین بار توسط 


۶ نشریه آب و خاک. جلد ۰۳۴ شماره ۲. خرداد - تیر ۱۳۹۹ 


مهندسین ارتش ایالت متحده برای مطالعات کنترل سیل حوضه 
رودخانه ماسکینگام در اوهایو توسعه داده شد. تانگ و همکاران 
(۱۹۸۵) براساس پژوهش‌های انجام شده. سه روش رگرسیون خطی 
(*1» شیب ترکیبی (20)) و داویدون - فلمنتچر- پاول (۳۳ظ) را بر 
پایه تحقیقات الگویی هوک و ژیور (111) برای تخمین پارامترهای 
مدل غیرخطی پیشنهاد کرده است (۲۴). موهان (۱۹۹۷) مدلی را 
برمبنای الگوریتم وراتتی (0۸) به‌منظور تخمین پارامترهای 
ماسکینگام غیرخطی ارائه نمود. نتایج نشان داد هیدروگراف جریان 
خروجی حاصل از روش 0۸ با هیدروگراف جریان خروجی مشاهداتی 
نسبت به روش‌های ارائه شده از سوی سایر محققان انطباق بالاتری 
دارد (۱۸). در پژوهشی بر اساس معادلات سنت ونانت روشی را برای 
روندیابی جریان غیر ماندگار ارائه نمودند که مشابه فرمول مورد 
استفاده در روش ماسکینگام با پارامترهای متغیر بوده و برای روندیابی 
از هیدروگراف اشل شده است. در این روش هیدروگراف دبی نیز به 
طور همزمان با هیدروگراف اشل روندیابی می‌شود و نتایج حاصل از 
کارایی مناسب نشان می‌دهد (۲۰). پرومال و همکاران (۲۰۰۱) با 
فرض ابت بودن شیب سطح آب در طول یک بازه کوچک از کانال و 
برقراری جریان ماندگار بین عمق در وسط بازه و دبی در مقطعی در 
پایین‌دست آن. روش ماسکینگام با پارامترهای متغیر را برای روندیابی 
موج سیل در کانال‌های منشوری با مقطع ثابت و جریان‌هایی که 
معادلات مقاومت بر آن‌ها حاکم است» بطور مستقیم از معادلات سنت 
ونانت به‌دست آورد. این روش در عمل از کارایی مناسبی برخوردار 
بوده و با حداقل اطلاعات نقشه‌برداری و بدون نیاز به واسنجی برای 
محاسبه زبری با تعیین پارامترهای روندیابی, نتایج قابل قبولی ارائه 
می‌دهد .)۱٩(‏ کیم و جیم (۲۰۰۱) الگوریتم ۲15 را در بهینه‌سازی 
پارامترهای ماسکینگام به‌کار بردند. نتایج به‌دست آمده از 115 تخمین 
بهتری نسبت به الگوریتم 0۸ داشت. روش آن‌ها نه تنها درکمینه 
کردن مجموع مربعات انحرافات (890) دبی روندیابی شده مشاهداتی 
و محاسباتی به‌عنوان تابع هدف موفق بوده» بلکه سایر پارامترهای در 
نظر گرفته شده مانند مجموع قدر مطلق انحرافات (8۸) دبی 
روندیابی شده مشاهداتی و محاسباتی و میزان انحراف دبی اوج 
مشاهداتی و محاسباتی (۳0) نیز مقادیر بهتری نسبت به سایر 
روش‌های پیشین بیان می‌نماید (۱۷). داس (۲۰۰۴) به‌منظور تخمین 
پارامترهای مدل ماسکینگام خطی و غیرخطی از یک الگوریتم تکرار 
شونده, بر مبنای حداقل کردن خطای هیدروگراف محاسباتی نسبت 
به هیدروگراف مشاهداتی بهره برد. نتایج نشان داد که بکارگیری 
الگوریتم. سرعت محاسبات را به نحو چشم‌گیری افزايش می‌دهد ٩[‏ 
چن و یانگ (۲۰۰۷) بهینه یابی پارامترهای ماسکینگام خطی را با 
استفاده از الگوریتمی مبتنی بر الگوریتم وراثتی را انجام دادند. کارایی 
این روش نسبت به سایر روش‌های روندیابی مورد بررسی قرار گرفت 
و نتایج نشان دهنده دقت بالای این روش را نشان داده است (۴). چو 


و چنگ (۲۰۰۹ الگوریتم 50 را برای تخمین پارامترهای مدل 
غیرخطی ماسکینگام با هدف کمینه کردن ۹50 به‌کار بردند. نتایج 
این مطالعه نشان از برتری ۳۹0۵ نسبت به 0۸ و اختلاف اندک 
نسبت به ۲15 داشته است [۸]. سامانی و شمسی‌پور (۲۰۱۰) برای 
روندیابی سیل در رودخانه‌های چند شاخه‌ای با بکارگیری روش 
پارامترهای مورد نیاز روندیابی از روش بهینه‌یابی غیرخطی (پاول) 
استفاده نمودند (۲۲). براتی و همکاران (۲۰۱۰) به ارزیابی عملکرد 
الگوریتم ۱۷1 به منظور برآورد پارامترهای دو مدل روندیابی غیرخطی 
ماسکینگام پرداختند. در این تحقیق پس از مقایسه عملکرد این 
الگوریتم ۷۷1 در برآورد پارامترهای مدل‌های روندیابی غیرخطی 
کارآمد می‌باشد (۱). محمدی قلعه‌نی و همکاران (۲۰۱۱) با استفاده از 
الگوریتم نورد شبیه‌سازی شده (9۸) دقت محاسبه تابع هدف 550 را 
با سایر انواع الگوریتم‌های مطرح شده تاکنون مقایسه نمودند. نتایج 
پژوهش آن‌ها نشان داد تخمین الگوریتم ٩۸‏ نسبت به تمام 
روش‌های گذشته به جز 115 از دقت بیشتری برخوردار است. همچنین 
این پژوهش گران در پژوهش دیگری در محاسبه تابع هدف 850 
نسبت به الگوریتم جستجوی هارمونی و فراکاوشی که در این زمینه 
به کار رفته است بهترین جواب به‌دست آمده است (۱۲). کارشان و 
رودخانه ویلسون و وای پرداختند و نتیجه گرفتند روش ترکیسی 
پیشنهادی در بین ۱۲ روش مختلف عملکرد بهتری داشته است (۱۵). 
ساقی و دلبری (۱۳۹۲)» پژوهشی بر روی عملکرد مدل‌های خطی و 
غیرخطی ماسکینگام جهت روندیابی متمرکز سیلاب با استفاده از داده 
های هیدروگراف ویلسون انجام دادند. در این پژوهش با استفاده از 
مدل 5۷/57۸7 و با بکارگیری اطلاعات مربوط به هیدروگراف‌هایی 
که با استفاده از نرم‌افزار ۸1۷18۳۸۴1 تولید شده بو پارامترهای 
مجهول مدل‌های خطی و غیر خطی به‌دست آمد. درنهایت معیارهای 
مختلف خطا برای مدل‌های مختلف محاسبه و نشان داد مدل 
غیرخطی نتایج بهتری به نسبت مدل خطی ارایه داده است (۲۱). 
حامدی و همکاران (۲۰۱۴) در پژوهش خود با بهبود مدل ماسکینگام 
غیرخطی از ترکیب دو رابطه ذخیره رایج استفاده کرده است. وی برای 
تخمین پارامترهای موجود در مدل ماسکینگام غیرخطی پیشنهادی از 
روش تعمیم یافته گرادیان مختصر شده (0150)) موجود در سلور 
اکسل و تابع مجموع مربعات انحرافات (۹90) دبی رونندیابی 
مشاهداتی و محاسباتی استفاده کرده است. بر این اساس استفاده از 
مدل ماسکینگام غیرخطی پیشنهادی سبب بهبود ۸۵ درصدی 550 
نسبت به سایر مدل‌های ماسکینگام غیرخطی در این مورد مطالعاتی 
شده است (۱۴). اکبری‌فرد و همکاران (۲۰۱۷) به بررسی پارامترهای 


هه یرای نا نادار اتکو ی رگن 
روی روخانه کارده پرداختند. آن‌ها برای بررسی عملکرد مناسب این 
الگوريتم از ششن تابع مک استاندارد با ایساد متفاوث با تحلیل 
حساسیت روی اندازه جمعیت مورد استفاده قرار دادند. نتایج حاصل از 
پژوهش آن‌ها نشان می‌دهد که الگوریتم مذکور نسبت به سایر 
الگوریتم‌های ژننیک» جستجوی هارمونی و ازدحام ذرات و رقابت 
استعماری برآورد مناسب‌تری از مقادیر بهینه پارامترهای مدل 
ماسکینگام غیرخطی نشان می‌دهد (۲). احترام و همکاران (۲۰۱۸) به 
بررسی معادله ماسکینگهام با استفاده از الگوریتم خفاش بهبود یافته 
برروی سه رودخانه ویلسون, کاراه ان میانمار پرداختند. این 
پژوهشگران برای بررسی عملکرد مناسب الگوریتم مذکور از هفت 
شاخص براساس فواصل تقسیم بندی شده برای افزايش دقت مدل‌ها 
مورد استفاده قرار دادند. آن‌ها نتیجه گرفتند که با استفاده از الگوریتم 
مذکور برای سیلاب‌های کاراهان ۸۱۱۷ بوده که در مقایسه با 
الگوریتم های مشابه نظیر ازدحام ذرات» ژنتیک» هارمونیک» خفاش به 
ترتیب ۰/۸۳ ۰/۸۸ 5۹۴ و ۱۲ کاهش یافته است. علاوه بر این, 
نتایج به‌دست آمده برای سیل میانمار نشان داد که ۹8 برای 
الگوریتم مذکور در مقایسه با ژنتیک» خفاش و ازدحام ذرات به ترتیب 
۲ ۱۱ و ۰۴۲ کاهش يافته و نتایج بهبود یافته‌تری را نشان 
می‌دهد (۱۰). بازرگران و نوروزی (۲۰۱۸) در پژوهش خود به بررسی 
معادله ماسکینگام با استفاده از الگوریتم بهینه‌سازی ازدحام ذرات برای 
رودخانه کارون در بالادست و پایین دست آن مطالعه کردند. آن‌ها 
مقدار دقیق خروجی جریان رودخانه را بررسی کرده و نسبت به رسم 
هیدروگراف خروجی اقدام نموده و نتیجه گرفتند اگر مقادیر ضرایب 
این الگوریتم متفاوت و غیر ثابت باشند دقت میزان دبی خروجی را 
افزايش داده و میزان خطای نسبی متوسط ۲/۴۴ بوده و در صورت 
ثابت بودن مقادیر ضرایب الگوریتم میزان خطای ۸۸۹ را ایجاد می‌کند 
(۲). فرزین و همکاران (۲۰۱۸) نسبت به روندیابی سیلاب رودخانه 
کاراهان و مقایسه آن با سیلاب رودخانه ویلسون و سیلاب رودخانه 
میانمار با استفاده از الگوربتم خفاش پرداختند. این پژوهش‌گران پس 
از محاسبات و تهیه هیدروگراف ورودی و خروجی سیلاب هر سه 
رودخانه و مقایسه الگوریتم مورد مطالعه خود با سایر الگوریتم‌ها از 
جمله ژننیک ازدحام ذرات. جستجوی هارونی میزان خطای نسبی 
برای رودخانه ویلسون به ترتیب ۸5 ۵۱ و ۶4 و همچنین برای 
سیلاب کاراهان به ترتیب ۸۸ ۸۳ و ٩۴‏ و برای رودخانه میانمار 
به ترتیب ۱۱ ۸۸ و 7/٩۴‏ کاهش يافته است و عملکرد این 
الگوریتم نسبت به سایر الگوربتم‌ها از دقت بیشتری برخوردار می‌باشد 
(۱۱). خلیفه و همکاران (۲۰۲۰) به بررسی بهینه‌سازی مسیریابی 
غیرخطی سیلاب رودخانه کارده بااستفاده از الگوریتم بهینه‌سازی ملخ 
(00۸) پرداختند و با سایر الگوریتم‌ها از جمله ژنتیک و هارمونیک 
مورد تجزیه و مقایسه قرار گرفت. تجزیه و تحلیل‌ها حاکی از بهینه 
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بودن اين الگوربتم نسبت به الگوریتم‌های مورد استفاده بوده و که 
مقادیر ۵/۲۹ ۲/۵۲ و ۵/۶۹ به ترتیب برای الگوریتم‌های بهینه‌سازی 
ملخ» ژنتیک» هارمونیک می‌باشد (۱۶). اکبری فرد و همکاران به 
بررسی پارامترهای معادله ماسکینگهام غیرخطی به وسیله الگوریتم 
جستجوی گرگ بر روی رودخانه کارده پرداختند. این محققین با 
مقایسه الگوریتم فوق و سایر الگوریتم‌های مورد استفاده در این زمینه 
از جمله الگوریتم گروه ذرات (۳۹0» جستجوی هارمونی (۳15) و 
رقابت استعماری (16۸) از شش تابع محک استاندارد با ابعاد متفاوت 
برای بررسی کارآمدی این الگوریتم‌ها مورد استفاده قرار داده و نتایج 
حاکی از آن بود که در میان الگوریتم‌های مذکور الگوریتم جستجوی 
گرگ فافزبه بآورد مناسب‌تری ازمق در رنه پاامتزضتین تال 
ماسکینگام غیرخطی اراثه نماید به طوریکه مقادیر ضریب تبیین و 
جذر ميانگین مربعات خطا برای رودخانه کارده به‌ترتیب برابر با 
۵ و ۰/۷۱۲۳۵۸ بدست آمد. 

باتوجه به بررسی منابع نشان می‌دهد که در سال‌های اخیر کاربرد 
روش‌های هوشمند و الگوربتم‌های فراکاوشی به‌منظور پیش‌بینی و 
بهینه‌سازی معادلات گوناگون مورد توجه محفقان علوم مختلف بوده و 
همواره نتایج مناسبی را به همراه داشته است. بررسی منابعی نشان 
می‌دهد که تاکنون الگوریتم جستجوی موجودات همزیست را محققان 
به‌صورت محدود به کار گرفته‌اند و کارایی این الگوریتم در زمینه 
روندیابی خروجی سیلاب تا اکنون ارزیابی نشده است. از این رو در 
این تحقیق به‌متظور تخمین پاامترهای مذل غیرخطی مانسکینگام از 
الگوریتم جستجوی موجودات همزیست 505 که بر پایه شبیه‌سازی 
رفتاری متقابل میان موجودات» که توسط ژنگ و همکاران درسال 
۴ ارائه شده است. استفاده شده است (۵). برای بررسی عملکرد 
الگوریتم توسعه داده شده نتایج حاصل از اجرای آن با روش‌های فرا 
ابتکاری دیگر مانند الگوریتم وراتتی (0۸ الگوریتم ازدحام ذرات 
(۴50) و الگوریتم رقابت استعماری (1۸) مورد مقایسه قرار گرفته 
است. کدنویسی الگوریتم‌های 05 ۸ ۳50 و 16۸ با استفاده از 
نرم‌افزار برنامه‌نویسی فلت نسخه ۲۰۱۸ نوشته شده است. 


مواد و روش‌ها 

مدل ماسکینگام از جمله روش‌های روندیابی هیدرولوژیکی است 
که اصول آن بر اصل پیوستگی جریان و معادله بین دبی و ذخیره 
موقت آب در طول مسیر استوار است (۱۲). اين مدل» روش به نسبت 
ساده‌ای بوده و دقت کافی در مسایل آبی دارد. این مدل به دلیل 
سادگی از بین انواع روش‌های هیدرولوژیک برای روندیابی سیل 
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کاربرد وسیعی دارد (۱۲). دو معادله پیوستگی و دخیره غیرخطی به 
عنوان معادلات پایه‌ای به صورت معادلات (۱) و (۲) در مدل 


1 

0 - [<ت- 

9 ۱) 
5 >] 2] +)1-26(6۵,[ ۲) 


در این معادلات ۰.٩,‏ ,1 و () به ترتیب ذخیره» ورودی و 
خروجی در زمان ) ۰ ک ضریب بدون بعد ذخیره برای رودخانه است 
که یک مقدار منطقی نزدیک به زمان گذر جریان از کل مسیر 
رودخانه می‌باشد (ثابت و بزرگتر از صفر)». 25 یک فاکتور وزنی که 
در مدل با روش ترسیمی به کمک سعی و خطا محاسبه می‌شود 
بین ٩,‏ و [,0(-1)+ ,2] هميشه خطی نیست. ویلسون و گیل 
در روندیابی سیلاب با استفاده از مدل غیرخطی ماسکینگام معادلات 
را به ترتیب مطابق معادله (۲) و (۴) به کار بردند (۱۳ و ۲۵). 
از 1-6(0۳)+ 201 ]1 ٩‏ 
(۴/ [1-260(0)+ ]22 5 

در مدل ماسکینگامپارامتر 0« به عنوان توان برای در نظر گرفتن 
اثرات غیرخطی (بزرگتر از یک) به معادله اضافه شده است. که مدل را 
قادر می‌سازد تا معادله غیرخطی بین ذخیره تجمعی و جریان را بهتر 
مدل‌سازی کند (۱۲). معادله (۲) نسبت به معادله (۴» به‌دلیل دقت 
بالات بیشتر رایج است و توسط موهان» کیم و جیم» چو و چنگ 
استفاده شده است (۸ ۱۷ و ۱۸). لذا با در نظر گرفتن معادله (۲) 
می‌توان معادله زیر را استخراج کرد: 


1 
1 و 9 1 
۵ ۱۱ مرت اک 
۵ و 0:2 
با ترکیب معادلات (۱) و (۵) معادله زیر بدست می‌آید: 
1 1 ۸ 
۶ 1 + ۳( - ۱ 
9 تا رح 
۸٩ ۰‏ ۱ 
که در آن سب تغییرات ذخیره نسبت به زمان می‌باشد. 
۸ 
۸٩5, ۳۲)‏ » 5 > ۳ 
که برابر ذخیره در زمان 1+ ات از آتجائی که با 


استفاده از روش سعی و خطا پیدا کردن مقادیر درستی برای ۰1 112 
و 26 دشوار و زمان بر می باشد» روشی متفاوت برای برآورد صحیح 
این پارامترها مورد نیاز می‌باشد. از این رو در اين تحقیق از الگوریتم 
نوین جستجوی موجودات همزیست استفاده شده تا این مشکا ر 


الگوریتم جستجوی موجودات همزیست ("808) 

الگوریتم‌های فرا ابتکاری جدید تقلیدی از پدیده‌های طبیعی 
هستند. الگوریتم 505 فعل و انفعالات همزیستی در رابطه بین دو 
گونه را شبیه‌سازی می‌کند. به‌نحوی که یک گونه به جستجو برای 
پیدا کردن مناسب‌ترین موجود می‌پردازد. الگوریتم 905 بایک 
جمعیت اولیه به نام اکوسیستم کار خود را آغاز می‌کند. در اکوسیستم 
اولیه» گروهی از موجودات (متغیر تصمیم) به‌طور تصادفی در فضای 
جستجو تولید می‌شوند. هر موجود زنده به‌عنوان یک کاندید از راه‌حل 
مسئله که با میزان برازش خاصی در ارتباط است» نشان‌دهنده درجه 
انطباق با هدف مورد نظر (مقدار تابع هدف) است. تقریباً تمام 
الگوریتم‌های فرا ابتکاری در هر تکرار» عملکرد جایگزینی را برای حل 
مسئله اعمال می‌کنند تا راه‌حل جدیدی برای تکرار بعدی به‌وجود 
آورند. در الگوریتم 5008 تولید راه حل جدید به وسیله تقلید تعامل 
بیولوژیک بین دو موجود در اکوسیستم اداره می‌شود. سه فازی 
همکاری (استفاده متقابل» هم‌سفرگی و انگلی که شبیه به مدل 
تعاملی بیولوژیکی در دنیای واقعی است معرفی می‌شوند. هویت هر 
تعامل مبتنی بر نوع تعامل» تعریف می‌شود. به این تریتب که سود دو 
طرفه بیانگر فازی همکاری» سود یک طرفه تصادفی با موجود 
دیگرتعامل می‌کند. اين فرآیند تا زمانی که معیار خاتمه فرآیند (رسیدن 
به حداکثر تعداد تکرار) برآورده شوند ادامه پیدا می‌کند. روند کلی این 
الگوريتم در شکل ۱ ارائه شده است (۵). 


یک مثال برای فاز همکاری که دو طرف از آن سود می‌برند 
یی که مرها تیوه سیگ ها ریا 
می‌کنند. شهد جمع می‌کنند تا به عسل تبدیل کنند؛ این کاری است 
که باعث سود بردن زنبورها می‌شود. این فعالیت همچنین باعث سود 
بردن گل‌ها می‌شود» چون زنبورها در طی روند کاری خود گرده را 
توزیع می‌کنند که گرده‌افشانی را تسهیل می‌کنند (۵). الگوریتم 505 
از این روند به عنوان فاز همکاری تقلید می‌کند. در "505 26 
همخوانی دارد. موجود دیگر که ,۸۶ باشد به‌صورت تصادفی انتخاب 
می‌شود تا با ,2۶ در تعامل باشد.هر دو موجود در رابطه تعامل با هدف 


افزایش مزایای بقاء دو طرفه در اکوسیستم شرکت می‌کنند. 


مدع فصعنصهعع0 ملامتصصررگ -1 


فاز همکاری 
ععهداین ]۳ 


2 


فار هم سفرگی 
روا وعصه‌حصحصوت 


ععصطحر 


قاز انگلی 
۱ 


شکل ۱ - روندنما الگوریتم 
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تولید جمعیت اولیه 
۷ص اه وال۳۵ 
صماهآن رم 


فرض بارامترهای اولیه 
۱۱۱۵( 


شرط توقف برقرار است؟ 


3 


[+ 


موجودات هم زیست 


۲ )۸۱0 505 -1 مساعز۲ 


راه‌حل‌های جدید برای ,25 و :26 بر اساس همکاری تعاملی 
۸ 
( ۷6۵۲۵۲۳1۲ -۱۷]۵۵۱۵۵1- 


2101)00,1(<>۴ ز 226 رریک2 


اوه 1۳0۳۷ 


( ۲۷6۵۱۵ میات ی ک1(<)2 ,۲۵0)00+ کر 

0) 

(۱۰) 
که در روابط فوق» (200)0,1: یک بردار از اعداد تصادفی است. 


2 


192" 


لت نامع ۵1-۷ ]۱۷ 


راظ ی 
" "و" به شرح زیر است. 

در طبیعت» بعضی از روابط تعاملی ممکن است سود بیشتری برای 
یک طرف از موجودات نسبت به طرف دیگر داشته باشد. به بیان 
دیگرء موجود ۸ ممکن است نسبت به موجود 13 سود بسیار بیشتری 
ببرد. در همین حال موجود ظ به سود خود بسنده 


می‌کند یا سود قابل توجهی نمی‌برد. در اینجاء فاکتور سود 


1 ۰ اب 
3 ] به‌صورت تصادفی در روابط (۱) و (۲) تعریف شده‌اند. این 


فاکتورها نشان‌دهنده سطح سود بردن هر موجود است که آیا یک 
موجود از اين رابطه سود کامل می‌برد يا بخشی از سود را می‌برد. 
رابطه (۱۰) برداری را نشان می‌دهد که به آن بردار متقابل گفته 


می‌شود و خصوصیات ارتباطی بین :* و ر > را نشان می‌دهد. قسمتی 
از رابطه ( ۷۵۵۲۵۲۲ - نت۱ ری 2) بو کاب سیری رنه 


به‌دست خواهد آمد. بر اساس فاز همکاری همه موجودات با افزایش 


اوعط 


توانایی و انطباق خود در محیطزیست روبرو هستند. بعضی از آن‌ها از 
رابطه همزیستی با موجود دیگری برای افزایش امتیاز بقاء استفاده 
می‌کنند. در اینجا هدف به‌دست آوردن م26 است» زیرا ,بمب دارای 
بیشترین درجه انطباق با نقطه بهینه سراسری می‌باشد. بنابراین» م2 
(راه حل سراسری) برای مدل کردن بالاترین درجه انطباق به‌عنوان 
نقطه نهایی و افزايش برازش هر دو موجود استفاده می‌شود. در نهایت 
موجود زمانی انطباق پیدا می‌کند که میزان برازندگی جدیدش نسبت 
به قبل از تعاملش بهتر باشد [۵) 
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فاز هم‌سفرگی 
رابطه هم‌سفرگی مانند رابطه بین ماهی چسبنده و کوسه است. 
ماهی چسبنده خود را به کوسه می‌چسباند و از غذای باقی‌مانده می- 
خورد بنابراین سود دریافت می‌کند. کوسه در برابر فعالیت ماهی 
چسبنده منفعتی نمی‌برد يا اگر هم ببرد به کمترین میزان است. مشابه 
با حالت تعامل دو طرفه, موجود :26 به‌صورت تصادفی از اکوسیستم 
مرتبط با موجود ,۸5 انتخاب می‌شود. در این شرایطء موجود ,کر 
حالی که موجود ,26 در این رابطه نه سودی می‌برد و نه متحمل 
ضرری می‌شود (۵). راهحل جدید ,۸5 بر طبق همزیستی هم‌سفرگی 
بین موجودات ,۸5 و 26 از رابطه (۱۱) محاسبه می‌شود. بر طبق 
قوانین حاکم بر طبیعت» موجود ,2 فقط وقتی انطباق پیدا می‌کند که 
میزان برازندگی جدیدش بهتر از برازندگی قبل از تعاملش باشد. 
قسمت ( رم میگ » به سودی که توسط :26 برای کمک به :26 به 
وجود می‌آید برمی‌گردد تا شانس نجات خود را در اکوسیستم کنونی 


7 ۸ 


امه 


,1-1 کر رک 


1060۷ 


فاز انگلی 

یک مثال از فاز انگلی» انگل پلاسمودیوم است که از رابطه‌اش» 
به وسیله پشه آنوفل از انسان به‌عنوان میزبان استفاده می‌کند. در 
حالی که انگل در بدن انسان رشند و تولید مفل می‌کنده فرد میزبان 
ممک ن است مالارب | گرفته و در نتیجه بمیرد (۲۵). 
در 05 موجود 26 به‌وسیله ایجاد انگلی مصنوعی به نام "بردار 
انگل "» نقشی شبیه به پشه آنوفل را ایفا می‌کند. بردار انگل در فضای 
جستجو با تکثیر موجود ,26 به وجود می‌آید. سپس تغییرات بر اساس 
انتخاب تصادفی و با استفاده از یک عدد تصادفی بر روی آن رخ 
می‌دهد. موجود :26 به عنوان میزبان و به صورت تصادفی از 
اکوسیستم انتخاب می‌شود. بردار انگل تلاش می‌کند تا 26 را در 
اکوسیستم جابجا کند. هر دو موجود ,»)با اندازه‌گیری میزان 
برازندگی (میزان تابع هدف) خود مورد ارزیابی قرار می‌گيرند. اگر بردار 
انگل برازندگی بیشتری داشته باشد. موجود ,> را خواهد کشت و 
جایگاه آن را در اکوسیستم تصاحب می‌کند اما اگر برازندگی 36 
بیشتر باشد» ,26 نسبت به انگل در امان خواهد بود و بردار انگل 
توانایی بیشتری برای زندگی در آن اکوسیستم را نداشسته و در آن 
اکوسیستم جایی ندارد (۵). 


مدل ریاضی 

در این پژوهش برای ارزیابی مقادیر بهینه پارامترهای که 26 و 10 
در مدل ماسکینگام تابع هدف به صورت کمینه‌سازی مجموع مربعات 
باقیمانده‌ها (850) بین حجم خروجی‌های واقعی و روندیابی شده 
مطابق معادله (۱۲) استفاده شده است (۱۵): 
۱/0۹۹0 
۳ 2 

که در آن ,06 حجم سیلاب خروجی روندیابی شده (محاسباتی) 
فرزنای اه راد کام‌های رسای روتفیانی فان نس اگیر 
مقادیر غیر ممکن 4 92۶ 11 انتخاب شوند. مقادیر 0( و 5 در 
مدل ماسکینگام منفی به‌دست می‌آیند و از اين رو از روش تابع غیر 
مستقیم جریمه به شکل زیر استفاده می شود: 


۱۲( 0 رربعک 1 تاره ات ر«بع5٩‏ 


0 «بع6 1 «اربی۵)|م0عوبی 6 
۱۳ 
7 0 2 بل . 2 ۳ 
که در آن و ثابت‌های جریمه بر و ببم6 مقادیر 


مثبت اما غیر واقعی و مقدار ذخیره و سیلاب خروجی در بازه زمانی 
1+ بعد از اعمال جریمه می‌باشند. 


معیارهای ارزیابی 

برای ارزیابی مدل‌های مورد نظر از شاخص کمینه‌سازی مجموع 
مربعات باقیمانده‌ها (850) و مربع میانگین خطای استاندارد (015۳) 
مطابق روابط (۱۲) و (۱۳) استفاده شد. 


کز0) ۲ 
9 سوورر 
۵ 0 

در روابط فوق ز5 سیلاب خروجی روندیابی شده, 5 ميانگین این 
داده‌ها (محاسباتی)» 0 سیلاب خروجی مشاهداتی. ۵ میانگین این 
داده‌ها (مشاهداتی) و 2 تعداد مشاهدات می‌باشد. ]1۷19 اختلاف 
بین داده‌های مشاهداتی و محاسباتی را نشان می‌دهد. که هر چه این 


مقادیر کمتر باشند نتیجه بهتری را نمایان می‌سازد. 


موّلفه‌های مورد استفاده 

برای ارزیابی الگوریتم های مورد بررسی در این پژوهش داده 
های رودخانه ویلسون(به عنوان مثال استاندارد) 9 رودخانه کارده 
(به‌عنوان مطالعه موردی واقعی) استفاده شده است. 


رودخانه ویلسون 
مثال مورد مطالعه در اين پژوهش برای اولین بار توسط ویلسون 
(۱۹۷۴) ارائه شده است. در اين منال رابطه غیرخطی بین مقادیر 


5 و[1>(0) + 2۲] توقرار قی‌باآشد وم قوان کرد 
الگوریتم‌های مختلف را در برآورد مقادیر بهینه پارامترهای معادله 
ماسکینگام غیرخطی مورد بررسی قرار داد. همچنین سایر روش‌های 
فراابتکاری در اين مثال مورد آزمون قرار گرفته است. بیشترین جریان 
ورودی و خروجی از اين رودخانه به ترتیب ۱۱۱ و ۸۵ متر مکعب بر 
ثانیه بوده است (۲ 


رودخانه کارده 

در این پژوهش برای بررسی عملکرد الگوریتم‌های مورد مطالعه 
در برآورد پارامترهای معادله ماسکینگام غیرخطی در مدل روندیابی 
سیلاب رودخانه کارده در استان خراسان رضوی به عنوان یک مثال 
واقعی مورد استفاده قرار گرفته است. حوضه‌ی آبریز کارده, با مساحتی 
معادل ٩۳/۲‏ کیلومتر مربع در شرق منطقه هزار مسجد کپهداغ واقع 
بوده و یکی از زیر حوضه‌های حوضه‌ی آبریز اصلی کشف‌رود می‌باشد 
(شکل ۲).بالاترین نقطهارتفاعی در شمال غرب حوضه ۲۹۳۷ متر از 
سطح دریا و پایین‌ترین نقطهٌ ارتفاعی در خروجی حوضه و در پایین 
دست آبادی کارده ۱۲۰۰ متر می‌باشد. متوسّط بارندگی در بخش 
جنوبی حوضه ۳۷۳/۲ میلی‌متر و در ارتفاعات بخش شمالی به ۴۵۰ 
میلی‌متر می‌رسد. داده‌های مورد بررسی در این پژوهش سیلاب دو 
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کاربست مقایسه‌ای الگوریتم جستجوی موجودات همزیست... ۳۷۱۲ 


روزه مشاهداتی در تاریخ ۱۳۷۱/۰۲/۰۷ تا ۱۳۷۱/۰۲/۰۸ به‌صورت ۶ 
ساعته می‌باشد که در دو ایستگاه کوشک آباد و کارده برداشت شده 


است. 


قبلا" بیان شد در مطالعه حاضر برای ارزیابی مقادیر بهینه 
پارامترهای > 26 و 0 در مدل ماسکینگام غیرخطیء کمینه‌سازی 
مجموع مربعات باقیمانده‌ها (550) بین حجم خروجی‌های واقعی و 
روندیابی شده به‌عنوان تابع هدف در نظر گرفته شد. نتایج پارامترهای 
بهینه به‌دست آمده از الگوریتم‌های مورد بررسی برای رودخانه 
ویلسون در جدول ۱ آورده شده است. در جدول ۳ بهترین مقادیر 
پارامترهای الگوریتم‌های مورد بررسی در مسئله روندیابی سیلاب با 
توجه به تحلیل حساسیت روی پارامترها نشان داده شده است. 
همانطور که در جدول ۱ مشاهده می‌شود الگوریتم‌های 505» ۳۹0 
و 10۸ توانسته‌اند 550 را به مقدار بهینه نزدیک‌تر کنند. شکل ۲ 
هیدروگراف خروجی روندیابی شده توسط الگوریتم‌های مختلف و 
هیدروگراف خروجی مشاهداتی برای رودخانه ویلسون را نشان 
می‌دهد. همانطور که در شکل ۲ مشاهده می‌شود تمامی الگوریتم‌های 
فرابتکاری توانسته‌اند به خوبی سیلاب خروجی را روندیابی کنند. 
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شکل ۲- موقعیت رودخانه کارده در حوضه کشف‌رود (شمال شرقی ایران) 
(صدو۱ 0۲ ماعم ۱۵) نموه 0مم/مصاعه معط رز رم زددام0 ها معط ۵۶ جمتای ما -2 مساهز۲ 


۲۳ نشریه آب و خاک. جلد ۳۴. شماره ۲. خرداد - تیر ۱۳۹۹ 


جدول ۴ مقادیر هیدروگراف خروجی روندیابی شده با استفاده 
الگوریتم‌های مختلف برای رودخانه ویلسون را نشان می‌دهد. جدول ۲ 
عملکرد الگوریتم‌های مورد بررسی در کمینه سازی مجموع مربعات 
باقیمانده‌ها (5500) را توسط پارامترهای آماری برای رودخانه ویلسون 
نشان می‌دهد. همانطور که در جدول ۲ مشاهده می‌شود الگوریتم 
5 در کمینه‌سازی مجموع مربعات باقیمان‌ده‌ها (650) با ۷51 
(۵/۸۵۳۹۲۹) و 550 (۱۲۸/۷۸۶۴۳۸) همانند الگوریتم‌های ۹0 و 
۸ عملکرد بهتری نسبت به الگوریتم 0۸ داشته‌اند. نتایج 
پارامترهای بهینه بدست آمده از الگوریتم‌های مورد بررسی برای 
رودخانه کارده در جدول ۵ آورده شده است. همانطور که در جدول ۵ 


مشاهده می‌شود الگوریتم‌های ۰505 ۳50 و 10۸ توانسته‌اند 550 
را به مقدار بهینه نزدیک‌تر کنند. شکل ۴ هیدروگراف خروجی 
روندیابی شده توسط الگوریتم‌های مختلف و هیدروگراف خروجی 
مشاهداتی برای رودخانه کارده را نان می‌دهد. همانطور که در شکل 
۴ مشاهده می‌شود تمامی الگوریتم‌های فراابتکاری توانسته‌اند به خوبی 
سیلاب خروجی را روندیابی نمایند. جدول ۷ مقادیر هیدروگراف 
خروجی روندیابی شده با استفاده الگوریتم‌های مختلف برای رودخانه 
کارده را نشان می‌دهد. 


پارامتر 


جدول ۱- مقادیر پارامترهای بدست آمده از الگوریتم‌های مورد بررسی برای رودخانه ویلسون 
۷۵۲ ۲۷۷۱۱۹0۱۱ ۲۵۲ عص)زدمعاج 0عوهصمه عمط رم 0عصنهاهان امه نو مط) 0۶ عمتاه۲ 1۳6 -1 12016 
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جدول ۲- شاخص‌های آماری مورد بررسی برای رودخانه ویلسون 
۵ ۷۵0۵ ع۱) ۲۵۲ عدم)هع01ه۱ ای‌نامتاداو -2 12016 


جدول ۲- مقادیر پارامترهای الگوریتم های مورد استفاده برای رودخانه کارده 
۲ ۱۱۵۸۲۰0۵۳ ۵۲ عصصط) تمعاه 0عومص0ه مط) رم 0مصتهاهاه )عم صصه ور مط) 0۶ کمتتاه۲۷ 1۳6 -3 12016 
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کاربست مقایسه‌ای الگوربتم جستجوی موجودات همزیست... 


جدول ۶- مقادیر روندیابی هیدروگراف سیلاب خروجی توسط الگوریتم‌های مدنظر برای رودخانه ویلسون 


۷۷۱۱۹0۵ ۲۳6 ۲۵۲ محصط)تدمواه معط وا 0عاهانهمای م2 0262 طامرمرهمنل چا 0مم۱ غناماجن معط ۵۶ معباه۷ 1۳6 -4 12016 
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دبی خروجی روندیابی نسده 
(متر مکعب بر ثانیه) 
(ورص) 001 0 عدناد1۳0 
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شکل ۳- هیدروگراف خروجی روندیابی شده توسط الگوریتم‌های مورد بررسی برای رودخانه ویلسون 
۵۲ ۲۷۷۱۱۹0۵۱ وه کم تمعاخ عم امرمهه0۵ ۴۱۲ ]م0 1 -3 ماعز۲ 


جدول ۵- مقادیر پارامترهای به‌دست آمده از الگوریتم‌های مورد بررسی برای رودخانه کارده 
۲ ۱۱۵۲۰0۵۳۴ ۵۲ فص تمعاه 0عومص0ه مط) رم 0مصتهاهاه ماع صصه ور م6 0۶ ومتزع۲۷ 1۳6 -5 12016 
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۴ نشربه آب و خاک. جلد ۳۴. شماره ۲. خرداد - تیر ۱۳۹۹ 
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شکل 4- هیدروگراف خروجی روندیابی شده توسط الگوریتم‌های مورد بررسی برای رودخانه کارده 
۲ ۱۱۵۲۰0۵۳۱ با عصصط)تتمعا۸ 60)م1 جامرحرهم0 و۲۱ )امن 1 -4 مساعز۲ 


جدول - شاخص‌های آماری مورد بررسی برای رودخانه کارده 
۵ ۱6) ۲0۲ م۲۱01 ای‌ناعتاهاد -6 12016 
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جدول ۷- مقادیر روندیابی هیدروگراف سیلاب خروجی توسط الگوریتم‌های مدنظر برای رودخانه کارده 
۷۵۲ ۱۵۲00 66 ۲۵۲ فصط)ز ماه عمط وا 0عافاناه مه 0242 مهم و 1000 )ناصایان م6 ۵۲ ومتااه۷ 16 -7 12016 


ی دبی ورودی دبی خروجی مشاهده‌ای دبی خروجی روندیابی شده (متر مکعب بر ثانیه) 
۹۹ "- (متر مکعب بر انیه) (متر مکعب بر ثانیه) (/2) 00 0 عصنات1۳0 
(ورس 1 و (وبیص 0۰۱ 0 ۵0 ۵۲و00 در ۵ ۳5 104 
0 2.41 2.41 241 241 241 241 
6 26 241 237 237 236 237 
12 37 241 131 156 121 136 
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جدول ۷ عملکرد الگوریتم‌های مورد بررسی در کمینه‌سازی 
مبوخ مرتمات تافیفانده‌ها (050) را توشظ پارامیرهای آماری رای 
رودخانه کارده نشان می‌دهد. همانطور که در جدول ۶ مشاهده 
می‌شود الگوریتم ۹05 در کمینه‌سازی مجموع مربعات باقیمان ده‌ها 
(550) با 25۴ (۰/۵۰۵) و ۹50 (۴/۵۵۲) همانند الگوریتم‌های 
0 و10۸ عملکرد بهتری نسبت به الگوریتم‌های 0۸ داشته‌اند. 
شکل ۵ میزان همگرایی الگوریتم‌های به‌کار رفته در رسیدن به مقدار 
مطلوب برای مسئله روندیابی سیل را نشان می دهد. همچنین این 
نشانگر همگرایی سریع ۹05 در مقایسه با سایر الگوریتم‌ها است. 

با توجه به برآوردپرامترهای غیرخطی معادله ماسکینگام توسط 
پژوهش‌گران در سال‌های اخیر برای رودخانه ویلسون نتایج آن در 


جدول ۸ با روش حاضر مقایسه و برتری روش حاضر بر سایر روش‌ها 


نتیجه گیری 

برآوردپاامترهای مدل ماسکینگام غیرخطی با استفاده از روش 
سعی و خطاء دشوار و با محاسبات طولانی همراه است. در سال‌های 
اخیر از روش‌های گوناگونی برای تخمین این پارامترها استفاده شسده 
است. روش‌های فراکاوشی از راه حل‌هایی بوده‌اند که توانسته‌اند در 
تخمین این پارامترها موفق عمل نمایند. 


کاربست مقایسه‌ای الگوریتم جستجوی موجودات همزیست... ‏ ۳۷۵ 


جدول ۸- مقایسه روش‌های مختلف محققین در سال‌های اخیر جهت بهینه کردن پپارامترهای غیرخطی رودخانه ویلسون 
که همع صز ممد ۲۵۵۵۵ ۵ کلمطامصه اصممنل ۵۶ مودصم -8 12016 
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در مطالعه حاضر الگوریتم نوین جستجوی موجودات 
همزیست(۹05) در تخمین پارامترهای مدل ماسکینگام غیرخطی 
برای سری داده‌های مشاهداتی دو مطالعه موردی شامل رودخانه 
عملکرد الگوریتم 05 نتایج حاصل از اجرای نش با الگوریتم‌های فرا 
گرفته است. در این پژوهش به‌منظور برآورد مقادیر بهینه پارامترهای 
مدل ماسکینگام غیرخطی در روندیابی سیلاب ابتدا مثال آزمایشی که 
به‌عنوان مورد مطالعاتی توسط پژوهشگران استفاده شده» تحلیل 
گردید. نتایج حاکی از بالا بودن دقت تخمین الگوریتم 505 نسبت به 
روش‌های گذشته بوده است. براساس نتایج به‌دست آمده از الگوریتم 


منابع 


مدنظر با استفاده از توابع هدف 5۹0 در مقایسه با دیگر الگوریتم‌های 
گزارش شده نسیت به۸» /۷ و نسبت به 10۸ ۰/۰۰۴ بهبود 
یافته است. در دومین مورد مطالعاتی که روندیابی یک سیلاب واقعی 
را نشان می دهد پارامترهای آماری به‌دست آمده برای رودخانه کارده 
توسط الگوریتم 505 برابر با 155 (۰/۵۰۵9) و 550 (۴/۵۵۲) 
بوده است که حاکی از عملکرد مناسب الگوریتم جستجوی موجودات 
همزیست در برآورد مقادیر بهینه پارامترهای مدل ماسکینگام 
غیرخطی در روندیابی سیلاب بوده است. در نهایت این پژوهش 
قابلیت الگوریتم 505 را در زمینه بهینه سازی پارامترهای مدل 
ماسکینگام نشان داد. از این رو استفاده از 05 جهت استفاده در 
زمینه روندیابی جریان با مدل ماسکینگام توصیه می‌شود. 


عناوم ۵۶ ممتمصصتاعظ مهد ۲0۲ صداتمواه وا ۵ ۵۶ امموع۲:۲ .2010 ۱۰ ۳۵0۵2۸۶2 20 ریک تاهتظ 220 ,.ن لته9ع( . -1 
صهتصه:1 -2010 .0۷ 9-11 ,رتمهم متایهیل و مهمتصم روممهعهگممن متانیهیل و۲ مفتصج طا9 ,۷۵۵8 مصیاممتاوت۱۷ 
(اعحعوهاه طفتاعص طزه صهتونهم م) بصعت رصحتطه 1 ,تفه نوا معین۷۵02 کمنهانه 1 .اع9۵0 ما۲ 

۲۱000 صام‌ملاوی ۲ تحمصنامم اه عطا ۶ه صمتامحطتاد۳۴ وتمامصصمت2ظ .2017 ۷۲ 20[مصصهتاه هه ,۷۲,۴ ۱۷۲۸020۲ ,رک تعللتهط اه -2 
یفصنهنا مه ممتاهع۱ .۳۷۵ طمتهک :تناو عععت) ره /۱۷) صطاتمعاه طمتقهعی ۷۷۵۱۶ وملهنا ۱۵061 عصتاتام۴ 
(امجمتواه طوتاعص طز صهتوع۲ 12) .95-112 :17067 موم مرمع متعم مامتا 

تقعصناً معط ۵۶ وعاعصصه‌تجم معط 10۶ فمبله۷ ماطمتنه۷ مصلوه 0۶ ]81160 معط ممتامع 1۳۷65۱ .2016 ۳۰ ۱۱0۲0۱21 مص2 ر.ل صقع227ظ . -3 
463-۰ :(32)14 )صمجصمم۵ه۷]2 وعمعتا0وم۴ عماج ۷۷ .(90ظ) تمه ندب مامتانهم مط عصتفنا 0مط)مصظ صیامصل و۱۷۲ 


9امعحطوطه طفتاعطط ره مدوم م) 


۶ نشریه آب و خاک. جلد ۰۳۴ شماره ۲. خرداد - تیر ۱۳۹۹ 


0۳2-0 هه 2960 1ع8مصظ صیاممف ۶0۲ ممتامصتادع معتعصصه‌تهم تقصصتام0۵ .2007 2 قصع۷ هه ر.ل صفط) 
849-۰ :(12)5 عمت)هاصنه. لهمتمصی مه ممممنهگ عدمصللصم صا فممتامم‌تص‌صصصمن حصط)تتمعله متاعصهع مصتاهتع 20081 
ماه . متاوتتمطمامصه .شمه خر بطمنوعد.. فصوتصهوعن . متامتطاصصنوو ...2014 .1۱۱ معمهعظ . هه .۱۲۷ عععط 
98-۰ :139 ماگ عک عمجم حصط) ت0۵ له 

کصمله تام مصونا معا هن ما عصتاندم۳۴ ۲۱۵۵۵ .2002 ,6.۷۲ عقکلانانهاظ مه ریکانک ۷۵ه]فه۷ نطو ربظ نامه 
413-۰ :(7)6 عصاتممص1مظ ملوم۲۵۱ ۱۷ ۵۶ امصینامژ ,«ماکصا صیتممافت ۱۷‏ 

۰ 12016101۱۲۶ ظ م1 ,۲۱۲02210116 عمجت موم0 .2009 ۷۰۲۰ «مط) 

عحعصنلصهماظ عطا ۵۶ ممتاحصتاوظ عاممصهتدوظ ما ممتاممتتاون) م9 مامتانظ عم رام۸۵ .2009 .۲:2 عصوطن 20 ر..۱ هل 
1024-7۰ :(14)9 مصیمعمتومط منوم۲۱۲0۲۵01 ۵۶ امصتم؟ ,۷۱۵061 یامملون۱۷۲ 

عصترممطصاعصظ ممهمتهنا لمح ممتامعن! ۵۶ مهم ,وام0ه۱۷]۵ یمن۱ 10۲ جمتامطصتاع تعامصه‌نجظ .2004 بر و12 
140-۰ :(13002 

.۸۵۱ عمط طوزه۱ .۲ ۱271( طفله۲ رب صوکه متعطامصلنل طخ وم ممعوع ۷ تمطل‌صت و۲ صمصطانم تقصاظ ری معط 
۶ ۷۵:۱ 2 عصلودا 00ظ)ععظ عصتاتا10 81000 حصیاع‌مناوت]۱۷ عطع عصدب0ع مرح .2018 بخ ملکمطفاط 4صه ,.۷ طمصنو.ط رطف 40تصفط 
۰ :(10)6 6۲اه 011)0.۷۷ع۵1 ۵2 4ص ممتاه2تصتامن هو ق2201 

۱ ۴۱2/0 هه ر.قل ۱۷۲۵۵0 .۷۲-۳۲ ۸۱۱271 ور تفلک بی صحعماداط رل تصقطهنه۲ وب تصهتفک .۷ دافم ریک ۲۵21 
۶ 100۲0۷۵4 مج اد 60منامي [86مجظر مصیتممیاوت تماممصهتهم عمط ۵ مصتوتا وعطموعع ونر صا ما۲۵ ۳1۵۵0 .2018 
0۰ :(10)9 2/6۲ 0210.۷۷ ع21 

۷۵۵۵ صاعمناوی ۱۷ عدعصتلمماظ عطا م۵ ممتامصصتاوظ ععاعمصهن۸ظ .2010 بک تصصتطه‌ظ مه رم ۲۵۸۵۵0 ع020۲ظ .۱۷ تصعلمطت 
(صهلوع۳۴ 12) .908-919 :(24)5 9011 20 ما2 ۷۷ ۵۶ لمصاوز ,عصتلهه‌ممم 2)60ا نله عمتوتا 

36)3-4(:353-۰ 020102۷ بط ۵۶ تمصریتان[.0مظ)مصه صنتعصتافت] ۱ عطع اطا عصتاتا10 ۳1۵0۵۵0 .1978 ,۸۷۲۰۸ 111 

سعص ح 0۷ عم صیممنادی۱۲ تقعمتلومه وبمتمحضا ‏ ,2014 هم طقط! صماع ۷ 20 رم ۳۸020 020۲8ظ .۳۲ ۲۱۸۲60 
0تطقطگ,2014 ۳۴۵ 29-1 )صعصمعمصها۱۷ مممسموم؟ عهاد ۷۷ چم ممصمنه‌کمم آقصمتنه 5۳ .صمتامناوه مهمتماه آمتممامصتطامی 
صحتوهم ع1) .صهیا - صهعطه ] . مصزممصتومظ مه ممصمزهگ دممعامومک عما۷۷/2 ۵۶ ممتاهژههوفوه صدتصه10-1و۵ نا تاطادعدهظ 
رامجمتدهاه مولعم از 

0 یام مامت تهممتمو ان مطع ۵۶ ممتاحفصتاع عامصهتع۳ .2013 7.۷۷۰ عع0) 200 و.ت) صقهآوتهتنات رب صقطه‌تفکز 
352-0۰ :(18)3 عمط منوما۲۱۲0۲0 ۵۶ امصمز معا تفع مرها 0زرها را 2 منوا [۱۷۲۵06 عصتاتام۴ 
تقمص جوا عمط ۵ ممتهمتصتانن) و0 قاعنا ,2020 بو هرطاخ 24 ویو فقصممطوهلمطک ویک تاتقصصوط ریق رفک 
30 مصتاا ۷ 161ظ ها 1۲2 گه تمصدمژ صطتمعلد همم عصلونا تعنک طمل‌ممک صا مادم ۲۱۵00 صصیامملوت۱۷۲ 
16/[.010.2020۰ 10.10 /1005://001.0۲2 

عصتعدا موه صصیامصناذنمهر تحعصتاممص فطا که ممتاحصتاوع ععامصصح‌تع۳ .2001 ,۴.8 حصتکا م20 .72.۷۷ موی یل نک 
131-۰ :(37)5 00120000وفظ وعمتبام وم م2 ۷۷ مهمتممصه عمط ۵۶ آقتامل 1۸۷۷۴۸ .1 اوعد رمق 

۶ مصدمل .صطاتتمعاه مناممه‌ی) عمتونا فلع۱۷]۵0 صیتومتافا تحمصتلمم ۵۶ ممتمصتافط تعاعصه‌ت2ظ ,1997 .و صقطم]۱۷ 
137-۰ :(123)2 عط1ممصلمظظ مل1ناهتل رز 

۶ 0۱2[ .۷۵۱۵0 حصیام‌صفن عامجصمیع۴ع۱ ۷2720 ۵ ۶ فومتاهء‌نان(۸ ۳۲۱۵10 .2001 .1:0 بازق۴ 20 ر.ظ۷۲ اقصییعظ 
196-7۰ :(603 عطهمصلمصظ ملو0۲۵۱۵ ۲۱۲ 

۰ 1 .۷۱۵۵4 متام طممهمل مها مامممنجع272۵۱ ۷ ,1998 .1,0 بازی؟ مه ری اقصیعظ 
115-۰ :(3)02 وطلهعصممط منعمامل ۲۱ گه امصتاه[ 

صامصلون ]۱۷ دبهم؟ عصتاندم: ۶۵۲ واممصظ تحمصنلمه مصه تقفصنا ممصمطمگهم ممتامع 18۷6۵ .2013 ۸ نداهن بط تطعمو 
طوناع‌ظ طانه صمتعع۳ ح) ,صقن ,عحواهه م82 :2013 ,66 17-19.واعنمم6 معتعوطظ صفتصمن؟ رعتصحصوظ فتبا۴ 15۳ ,08م1] 
(2050201 

تحعصنلطمه هز فصصماولاه 1۱۷۵۲ 0عطمصنهعه ‏ ظا عصتانام: ۶1۵00 ۳۱۷0۲۵10210 .2004 .0 عناممتقصفطاه. مضه ۳۰۸ تصمهو 
55-۰ :(42)1 ظ6270وم۴ متلییهیل ۲ گه تمصیامول بممتاهمتصتامه 

ملظ ۶ه تمصم! .عصقنم: ۶1۵00 موی تهممتامومه گ۵ فلورلممه ,1987 ۲.۲۱۰ وماهاندعد 4ص ر.ظ.۷ طومنه 
2-۰ :(۱13)1 ۲166۲1 

:(111)12 مدزمعصزمص۴ متلنییل ۲ ۵۶ امصصاول .0مصامجظ مصیام‌متافی ۱۱ عحعصتلممه ها عصتانام۲ ۲1000 عبات :1985 بک۰ ۷ عصتا 1 
1447-0۰ 

۵ ۱۲۵۱۱۵0۵ رمتتطاوم‌صها! رومتاهعن۲0 مه۱۷۱2۵۱۱ ,وومامتهبظ ممتمعصتمو۲ 1974 ۲۰۷۲۰ طموا1 ۱۷۷ 

حمن ناوت ۳۷2۵۲ ۵۶ ط0ه2تصتان 0۵ :2019 ۷۲۰ تصهصطرمک عقصصمصطمم مظعم و ۷ ام11هطک ریم تصعوظ رو 121116 
۶ فتاه .۵0۷02 ۵ عه۵ :9000۷۰ 286 081مه قمطع2۸۲۵۲0 ۷۷ طا فصطتمعله ملامصمع عصلوتا فصظماو 
(اعحعوط9اه طوذاع ۴ ۷16۲ صهزورع۳ م]) .111-124 :(32)3 ۱۵۸۵186۵۵1182 01۷11 ۲۵۲00۱۷۷5۹1 


25- 
26- 


نشریه آب و خاک (علوم و صنایع کشاورزی) م6 0ص ۲۷26۵ ۵۶ تقصتته[ 
جلد ۰۳۴ شماره ۲. خرداد - تیر 1۳۹۹ ص. ۳۶۵-۳۷۸ رد 365-8 ۲۰ ,2020 6 [- 1۷۲2 ,2 ,0 ,34 ,۷۵۱ 


]۱۷۱60۵۲ 6 )معا طمسععی معتصمعه۵ »تمتصاصصوی مه ۵۶ ممعتنجمرهمن 
۷۲00۵ 0۱۵۱۸۵ ۲۱۵00 دا محصطااتتنمعاخ 


"عقصم‌طوه م1 .6 - تاتفصصوظ :16 اتقو هیک - طمگن(فط)1 .6 
05-02-2020 :16661۷60 
18-03-2020 :6660۲60 


,۷ طا دعتالمناعمه 1 صرح لمتمصمصت منامتمصصناهر وفنیقع مه عغقطا مممم‌جورمصمطم لقسقهه ه ود ۳۱۵۵۵ :ممتام00 1۳۸۲۵ 
ص 20 ۵۵00 ۵ عصتسل لقع ولا مه میا ومتافتهام هه و1 عص ۱۵ رمعمگمعمط 1 .قهمیه لهتتانهمتتيه مضه رلفتتاون0ط1 
وه هه ومتافتمامهتهصل ففقطا هصتامهعع10۲ رلحعصهع ما .عمصه)ع0ممح1 اومصصانا فص ۵۶ ود ۲۱۷۵ عمط ۵۶ وعمهاض اصعع01 
۶۲ 3 0۶ م0 و1 عص1۳۵ ۳۲۱۵۵۵ .05مطامم مادام 11000 مطا عصتونه ففمل 9۵ 10اهم ۲۷۵۲ مطا صر تما 0صه طامع0 10 
0 وه1۱۷ 0۳21060 1 «اتهلنمتاتهن ممصصصرمی ود عصتانام؟ ملعمام۳۱۷۱ مامعز۳0 عصمعصتهصه م۳2 م1 ومتاوع1 ]1۳0001 
ما ممتصصعاعل مه مممتامعو-وو0۵ م0۵ 10 8660 2 15 06268 رتقطا مل 10 باق صلفهه-۱0 0۶ عامجا عمط ۱1 ۲۱۷۵۲5 
۱۱07۲0102۱ عطع صر ولمم منام۲ مه 0عویا «1عصمناوع ور 0مطامصه لاف م1 بمتوع۲ وه ما 5100۵ ۲1۷۵۲ 
,۲۳6۵۲6۵]016 ۲۵۵۵۲0۰ وتا صا ملناوع پرتمامه]فتاهه م«مداو عقط وعمطاممور مصط)تمعاه متافتتنمطهاممه ۵ ممتاهاجممصماموصظا مط ]1 
لمصتاوه عمط عصتامستاوه ما (رقمق) همع فحصوتصهع:0 متاماصصیود مطا که پمصمته‌تله مطا 0مامناهبه ۵ ,تاه فلطا ظ1 
01۰ نا ء مافت ۷ تهممت- عم عمط ۵۶ ممتامصتاوع تمامصصهته0 

(د0ق) همع مصعتصهعنن0 متامتطارد ۵ ممصهصمهن فطا 0ماهتلهه بای ونظ 1 ءفل‌مط)ع۷۱ مضه واهته) ۱۷2 
معط ماهعناوع1۱۵۷6 10 ,۵001 تهممت مها صیاوملافی معط ۵۶ فعمامصونهم صتصتاوه فطا عصتاقصتهه ما صصاتتمعا2 
فص ما مومع م۵ (وم0ق) تفع فصصعتصعع0 متامتماصصروی ما ۵۶ فالناوع مطا وطتلاطاه‌تنوع وعصتصق و طم‌تمعوم۳ 
016۲ 01 عاانت۱ 
معط 20 ,(ر0وظ) ممهمتتانن0۵ ج ۳ م۲2 معط ,ررهم)) صط)تعمعله عتامصمی فطع عصتلتامض1 م0مطاممه متافتتناع]-ها۷]6( 
6۵ 10 12726 1۵0 عع2 طمنط/< فصمتانتامو ۵۶ اعد 2 و1 مهد ومتادزم-۱۷۲۵)۵ .ره) صطاتتمعاه 0۷۵ا۵ص۵) )و1ل۳0۵۵۲12] 
۹01۷60 عصله0 مامت ممتاهم2تصتاونه مطا غبامهاه فصمتام‌حصاووج 1۵ ماقم رفح دمتاوتتمصماه]۱ .160مصصهو «اعاعآم‌ممع 
۵ مطصمنامهع)ص1 مط) وماهاته صطاتتم‌عاه قو .عصمامامعم ۵۶ بامتنه۷ 2 10۶ عاونا 0۵ ,876010۲6 ,مج رمطا هه 
لمتانصاً صه طاذ عانتما صطاتموله قمو .ماطد‌اننه )ومع فطع فص ما اوعد فعتهعمو عصم قطا ۷۵۲ 2 طاً ومل۹۵۵۵ 0۷0 
مه (ماصمتیه۲ ممنوتهع0) مصصوتصهعته 0۴ نامع 2 رد600 ]0 فمعهاو اوه مط 1 دومع 21160 ممت)هل0۲ا20۵ 
صتقاهی 2 10 ۵0905وع۳هه فص ممنتا۵و 2 1۵۲ 10246صهع 2 و1 مصوتصهوه مق .عمهمو طمتوعو مص ما 0م2تمصمم «امهمه1۵ 
00۰ص 00600۷۵ 0۶ اصیامصصع) عوممرتام 060ص مطا ط نونجم ۵1 معتوعل مطا عصتامهمع:۱6۵ بات 0۶ 062۲62 
ع ص وملهعوو 0۵0 مطا ممعهه طمتامهعاص لهمتومامزه معط عصتام‌تصتص رها متام فص ه قمویا صط)تدمعله وتط ]1 
۵ 010102102 عطا م6 تقاتصته رمتاتوهته 24 رصوتاه‌عممصصصم ررمف-وومع) وعفقطام )مصتافنل ععتط1 ,دودعم 
,۷۷۵۷ عتطا ما .صمتامحتماصا ۵۶ ۷۵۵ مطا وم 0عووها 0160 و1 متامجماص1 مق ,م1000 و27 ,۷۷۵۲1 ۲۵21 6ظ6 صرّ ۸۵06 
40 ,00256 صوللهممم‌صصصمی مطا ماصموم۲م۲۵ ۵۲01 7۵۲ م0 مط) رفقمام ممتهعو6۵0 فطا فامموم۲۵۵1 ]۵۲0 ۷۵-۵۷ عط 
0 0۵1۸8۵6 15 موه روعفقطام 21 م۲ .عفقطم متالفهتهم معط کصعوع۲۵۵۲ وعووم1 0و تعطان فطه فصه مه زهسمون فطا 
۲اه نحل تقنظ ما عصتطمهع) 0عامامصمم و1 ووععمم فطع مایا ومیمتاومم قوععمقم فظ1۳. ,تمطان فطا 6 «اصهم‌صه1 
0 معصهاوض1 ۲621 2 فه معومهام ۷۵۵ 0۲0۷1۵۵۵ 1۵22۷1 صه0128ظ صد بع ۱ طاملعقکا فطع رطم‌تهعوع فطع ما .(فصمتاه1۵ 0 
و ۱۱ 
۵ (890۵) هنانز ۵۶ حصتی فطع مط21تصتصتمظ رای فتط ق .1ع0مصظ هصتنامه ۶1000 عط صا ممتامناوه صامماوی۱۷ 
تمه مطا ماقفبلهه مه ممتاعصیظ متاععژه0اه صج فج 0۲۵0م0اقصمي فه نامام ما۲0 20 ۲۵2 ۵۶ مصصتتان۷ معط 
همهم اون 0مصنهاه9اه م1 مصمتفنوه صتعصتمی تهمصتومه فص هد قظر مضه کم یک ۵۶ 0272۳66۲ 
لقصتامه ما افو معط مافصصتدم؟مم2 لیام فصصط)تتمعاه 10۸ 20 ,ر0وظ روم مطا عقط) 0760و عبر طامها 10 وحصط) تع0ع21 
۰ 2۶5 11000 )ناوتان معط ما۲0 لیام فحصط)تزمعله متاوتتام-هاممه 21 20 میتله۷ 

۶ ده مطا ۵۲ صمتاه2تصتمتمه ما ۵0مطو هنز معط عم قالناعع؟ صطاتنمعاه مقو م1 خصمتععیه عزظ 20 واناعع1 
۵ ۱۷۲9۳ 20 رتع۷ ۳۲ ۱۷۷1190 عمط ۵۲ 128.78 ۵ لقناوم وه مود 280 5.85 ها لقناوع فهب او مت (ریقوه) وعتقتا0٩‏ 
0 مط) مه محصصهه مصطاتتمعله هت عمط صقطا باتلمصمتامصنط تمتاهه 24ظ هه 4.55 ۵ لماوع 7728 فقو 20 0.505 م6 لقناوم 
معط متقصتاده 10 ۵060ععممو طمنطه فصمتانامو ما ت۷۵ 09مطامص. متافزاع-مامظ م1 . .فصطاتتموله م۲۸ هه 
عومص-ع۵ظ فطا متقصتایی مه 9۵0 ۵۵ (قم0ق) طمنععه فصوتصهع:ن متامصاصررو امه فطع ,تفای فنطا ق .عتمامصهتهم 
۵ ۳/۷۵۲۵ ۷۷۱۱۹0۵ هه م۵0 ۵ ومتهناه هنز وه 0 ماج تقممتاه وان م1 ,قتماممصهتهم ام0مصظ مصتاعملآفت] 


۳۲۳۵1۵560۲ مرح 0۲ووع]۲۵ظ 001246وفظ مناد ۷۷/2۲ اصم0ناو مق رمتطمدااگ اج ۱۷۷ 010216صون) ,۳۵,۲۸ 4 0صع 3 ,2 و1 
۱ 

(1. ۱.26 ۵ 92۲11ق: اتمصظ میمصت طمم دمم -۶) 

1701: 10.22067/157.۷73412.85190 


۸ نشریه آب و خاک. جلد ۰۳۴ شماره ۲. خرداد - تیر ۱۳۹۹ 


0۸۵۱ که ماو فحصطاتتمعله متاوتتتعطمامص تعطانم مه 0عتومطصمن ۸۵۲۵ صمتاه]ممجصمام‌مصا ق0و ۵۶ فعاناوع: م1 .۵0«مامرصرم 
۱ ۱۵۵60 عامرجصقه آم)ممجتنومه مطا رولاعن رجم‌نهعوعر ونط ما .وانلمدمتامعصیظ دمو مطا مف2عتاوع1۳۷7 10 ۲0۸ 4صه ,0وظ 
۰ ۴۱0061 ت۲0 11000 صعمتام ]۱۷ توممتا-مماظ تاو ما واه م6 0مهع10۷6901 ۷۸۵۵ وتعطمتهعوع1 فص 
عص 0 عصت:0ععظ .08مصامصظ فنامدعيم فطا ما ماطمهم‌مه مه عمتتقصتاده وفامعتم 90۵ فمط 060و فالنافع: مط 1" 
2 4 محصطاتمعاه ۲۵۵0۲۲۵۵ تعطان مه مه 0متهم‌جصرمع قه فصمتامصظ ۵اععزهاه فط) عصئوت وه ۲0۷60مصص1 قهع/ م9 مطا رفلتاوع۲ 
اوه رعصتاا10 11000 ۲۵21 2 ود فنص رتم۱۱۷ آمامم‌صصتهمره 9۵6080 فص صز مه ط 0.004۵ 220 یه ط1 796 0۴ ۲۵/0 
کنط .صطاتمواه دو معط مر موق 101 4.5528 20 ۱۷۲5 ۶0۶ 0.5059 ۱۷۵۲۵ ۳۲۷۵۲ امعم عمط ۲0 تماعصصه‌تدم ادمتادتاماه 
صنصتامن فطع مصتامصتاوع صا مصطاتتمواه طمعمعو محصعتصهعر0 متامتصصصررگ مط ۵۶ هتکمممتامصنظ مامتع200100 ما عقظ) فمطه 
0 «انلزاه و معط فاطمتاطعنط ۵ طمهعوع: متا ,لفط .عتمامصصه‌تهم اع0مصه مصتاتا0؟ 1000 مصیام‌متاوی ۱۷ تحعمتا-ممل 
۰ 0161 صیعمیآوت ]۱۷ عط) ممتصتامه 

۰ 10۸ 20 0وظ مط قه میاه رم‌صممموج 00مع 0760و ۵و رصهع‌تاو 11000 م99 معط مد خصمتعسمومن 
م6 1۳0۲0۷۵0 هه ,اف هه موه صمتاعصظ ع۲تامعزد0اه مطا ص هت فمطا صقطا ماع 130 و۷2 ۵و بلعهعع1 ونطا 12 
عط معامجع0 فنط1 ,۲0۸ 0صج 0وظ معط ما آعع۵و۲۵ ظا ,۲6۹0۵۵۵۷۵ ,466 200 0.002 و۵ ۱۷۲9 4ص فقو صمتاعصنظ 00600۷۵ 
تحعص-صمه حصاحصتامن فطع عصتامحصتاوه صا صاتتمماه بقع فصعتصمعنن متامتطصروگ مطا ۵ رتلمممتامصنظ ماقت۲ع0۲0م2 
ناموت۱۷۲ مطا مم2نصتامه ۵ اتلتداح دم مطا ععمعت4ص فعصتع‌ص م1 ,قمامصصه‌تهم [06مع مصتادام۲ 1000 ناموت۱۷۲ 
1۵۵0۵۰ 18 [6مظ مام‌صیآدی ۱۷۲ قط ط16 عصتاتا10 ۶1000 طر دوه معط مصتوتا رععمگع۵۲ظ۰۲۳ .وععامصصه‌تهم 2۵061 


,۵06 صصیام‌صمن . ,متام ه)ه۱۵ رمتافتتمطفاها متام ملومل۵0ع0 ۲۱ روته ۱۷‏ 121060 :1۵۲۱۷۵۲9 
هام0۵ 


